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ABSTRACT 

We report on spectroscopic observations covering most of the 475 BL Lacs in the 2 nd Fermi LAT 
catalog of AGN. Including archival measurements (correcting several erroneous literature values) we 
now have spectroscopic redshifts for 44% of the BL Lacs. We establish firm lower redshift limits 
via intervening absorption systems and statistical lower limits via searches for host galaxies for an 
additional 51% of the sample leaving only 5% of the BL Lacs unconstrained. The new redshifts raise 
the median spectroscopic z from 0.23 to 0.33 and include redshifts as large as z = 2.471. Spectroscopic 
redshift minima from intervening absorbers have z = 0.70, showing a substantial fraction at large z and 
arguing against strong negative evolution. We find that detected BL Lac hosts are bright ellipticals 
with black hole masses M. ~ 10 8,5-9 , substantially larger than the mean of optical AGN and LAT 
Flat Spectrum Radio Quasar samples. A slow increase in M. with z may be due to selection bias. We 
find that the power-law dominance of the optical spectrum extends to extreme values, but this does 
not strongly correlate with the 7-ray properties, suggesting that strong beaming is the primary cause 
of the range in continuum dominance. 

Subject headings: BL Lacertae objects: general — galaxies: active — Gamma rays: galaxies — 
quasars: general — surveys 



1. INTRODUCTION 



The Fermi Second Source Catalog (jNolan et al.ll2012l 
2FGL) lists the 1873 most si gnificant sources de tected 
by the Large Area Telescope (jAtwood et al .1120091 LAT) 
during Fermi's first two years of sky survey observations. 
The majority of these sources are associated with jet- 
dominated Active Galactic Nuclei, the so-called blazars, 
many of which are bright, compact radio sources. There 
are, in fact, 1121 such associations (1017 at \b\ > 10°), 
collected in the Second Catalo g of AGN Detected by the 
Fermi LAT (jAckermann et al.ll2QllL 2LAC). These AGN 
are further classified as Flat-Spectrum Radio Quasars 
(FSRQ) where the optical spectrum is dominated by 
thermal disk and broad-line region emission, BL Lacs 
(BLL), where the optical spectrum is dominated by con- 
tinuum synchrotron radiation, and a collection of miscel- 
laneous, mostly low luminosity related sources. In 2LAC, 
the sample included 410 BLL, 357 FSRQ, 28 AGN of 
other type (generally low z, lower luminosity Seyferts), 
and 326 AGN of (then) unknown type. 

These 'Blazars' (BLL and FSRQ) are the brightest 
extra- Galactic point sources in the microwave and 7-ray 
bands; study of their population and evolution are cen- 
tral topics in high energy astrophysics. To support such 
studies we have acquired sensitive spectrosco pic observa- 
tions of this sample. In a companion paper (jShaw et~aT1 
I2012L hereafter SI 2), we reported on measurements of 
a large fraction of the FSRQ. Here we concentrate on 
the BL Lac objects. Our study has also found types 
for some of the unclassified blazars; the 'unknowns' have 
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now decreased to 215 (19%), and the confirmed BLLs 
have increased to 475 (42% of all 2LAC AGN). 

In $21 we outline the sample properties, data collection, 
and data reduction steps. We also summarize principal 
features of the spectra. In $3l we describe our spectro- 
scopic constraints on the redshift, including a technique 
to provide uniform redshift limits based on searches for 
host galaxy emission. In $11 we give estimates of the BLL 
black hole masses. We turn to comments on the principal 
BLL feature, the non-thermal dominance in the optical 
in $5j and conclude with general remarks in $6l 

In this paper, we assume an approximate concordance 
cosmology - Vt rn = 0.3, = 0.7, and Ho = 70 km s _1 
Mpc -1 . 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. The BLL Sample 

BLLs were originally identified as optical violently vari- 
able AGN, and are often characterized by an optical 
continuum dominated by synchrotron emission. Their 
broad-band spectral energy distribution (SED) is de- 
scribed by a synchrotron component peaking in the far- 
IR to X-ray bands and an Inverse Compton (IC) com- 
ponent peaking in the MeV to TeV range. In the radio 
these sources display str ong core dominance. A ccord- 
ing to the unified model (jUrrv fc Padovanilll995[ ) BLLs 
are the beamed counterparts of the FR I radio galaxy 
population, while the FSRQ are associated with FR lis. 
However, the principal BLL characteristic, a dominant 
and varying synchrotron/IC continuum, is a sign of a 
powerful jet whose emission is beamed closely toward 
the Earth line of sight. Thus the distinction between the 
traditional BLL and the FSRQ is se nsitive to the precise 
state and orientation of the jet (e.g. iGiommi "eTaLll20T2h 
and, indeed, variations in jet power or direction can bring 
individual sources in or out of the BLL class (S12). 
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Our own assignment of the BLL label f ollows a prag- 
matic "optical spectroscopic" definition (Mar cha et al.l 
1996): these are blazars with no emission lines greater 
than 5 A observed equivalent width, and a limit on any 
possible 4000 A spectral break of < 40% ([Marcha et al.l 
119961 : iHealev et all 120081 : IShaw et al.l 120091 ). For transi- 
tion objects with varying continuum emission, we retain 
the BLL label if it has ever been confirmed to be a BLL 
with a high quality spectrum, even if subsequently ob- 
served in an FSRQ state. Within the BLL population 
it is common to classify sources based on the peak fre- 
quency of their SED's synchrotron component, as esti- 
mated from radio/optical/X-ray flux ratios, separating 
the sources into low-peak {v pea k < 10 14 Hz, LBL), in- 
termediate peak (10 14 Hz< v pea k < 10 15 Hz IBL), and 
high-peak (v pe ak > 10 15 Hz, HBL) sources. We adopt 
here the LBL/IBL/HBL designations from 2LAC, which 
provides such subclasses for 74% of all 2LAC BLL, and 
83% of those with spectroscopic redshifts (Table [T]). 

The evolution of BLL has long been controversial, and 
it has been claimed that they are predominantly a low 
redsh ift population, showin g strong 'negative' evolution 
e.g. iBeckmann et al.ll2003h . especially for the HBL class 
Rector et al.ll2000f ). One challenge to studies of the cos- 
mological evolution of these sources is the difficulty of ob- 
taining redshifts from their nearly featureless, continuum 
dominated spectra. Indeed, many of the early studies 
using X-ray or radio-selected samples had highly incom- 
plete redshift measurements, even though the samples 
were confined to relatively bright sources. Uncertainty in 
extrapolating from the measured set of redshifts compli- 
cated population interpretations. In the Fermi era, this 
issue becomes critical, as the large BLL contribution to 
the LAT source population and the hard BLL 7-ray spec- 
tra ensure that these sources are a major fraction of the 
cosmic 7-ray background and may, indeed dominate the 
LAT background at high energies (Ajello et al, in prep). 

Since the LAT provides a large, uniform, sensitivity 
limited (~ flux limited) blazar sample, it provides a new 
opportunity to make progress on these questions. Im- 
portant to interpreting the LAT blazars are the strong 
correlations between the LAT-detected (IC) part of the 
SED and the synchrotron component covering the op- 
tical band. The synchrotron peak location determines 
the sub-classification, but also correlates with the inten- 
sity and LAT-band spectral index of the IC component, 
which in turn affects the depth of the LAT sample. Fur- 
ther, as illustrated in this paper, the synchrotron peak 
and intensity also affect the difficulty of optical spec- 
troscopic measurements. Since we wish to recover the 
detailed evolution of the blazars, preferably also follow- 
ing differences among the LBL/IBL/HBL subclasses (or 
even better, a physical parent property leading to these 
subclass es), one need s a detailed treatment of both the 7- 
ray (e.g. lAjello et al.ll2012af ) and optical selection effects. 
We reserve such analysis for a future study, noting here 
only the most prominent trends in the measured optical 
sample. Of course, characterization and minimization 
of the optical biases are greatly aided by high redshift 
completeness, the goal of the present paper. 

We have accordingly studied the 2LAC sample over a 
number of years with a wide variety of telescopes, striving 
to be as complete as possible. We here report on obser- 



Table 1 

Sample Completeness 



Set 


Total 


Spec z 


z min Unk. Z 


2LAC 


1121 






BLL 


475 


209 


241 25 


LBL 


72 (21%) 


35 (20%) 


36 (23%) 


IBL 


91 (26%) 


41 (24%) 


49 (31%) 


HBL 


187 (53%) 


98 (56%) 


73 (46%) 



Note. — BLL includes 65 AGN so classified since 
the 2LAC paper. 349 BLL have 2LAC SED sub- 
classes; percentages give the breakdown. 174 of the 
spectroscopic redshifts and 158 of the lower limits 
have subclasses; percentage breakdowns are given. 



vations of 278 BLL objects and 19 other Fermi blazars, 
not included in S12. We further analyzed 75 SDSS spec- 
tra, treating them in the same manner as our new ob- 
servations. Note that, unless we suspected a published 
redshift was erroneous, we generally did not obtain new 
spectra of many of the brighter, famous BLL with red- 
shifts in the literature. In several cases when we obtained 
new data it strongly contradicted the literature redshift, 
either from a new secure spectroscopic z or from an inter- 
vening absorber at larger z. In the end we retained 107 
redshifts from literature values; however for several early 
BLL z's we have only inspected plotted spectra (of vary- 
ing qualities); we suspect that at least a few erroneous 
values remain in this set. Our new data provide new 
spectroscopic redshifts for 102 objects and secure lower 
limits for many more, as summarized in Table [TJ This 
brings the 2LAC BLL sample to 44% redshift complete- 
ness, and 95% completeness including redshift limits. 

2.2. Observations 

The quest for high completeness has driven us to em- 
ploy medium and large telescopes in both hemispheres. 
Observations were obtained with the Marcario Low Res- 
olution Spectrograph (LRS) on the Hobby-Eberly Tele- 
scope (HET), with the ESQ Faint Object Spectrograph 
and Camera (jBuzzoni et all [l98i EF OSC2) and ESO 
Multi-Mode Instrument (jDekker et al.lll986l EMMI) on 
the New Technology Telescope at La Silla Observa- 
tory (NTT), with the Goodman High Throughput Spec- 
trograph (GHTS) on the Southern Astrophysical Re- 
search (SOAR) Telescope, with the Double Spectrograph 
(DBSP) on the 200" Hale Telescope at Mt. Palomar, 
with the FOcal Reducer a nd low dispersion Spectrograph 
([Appenzeller et al.l I1998L FORS2) on the Very Large 
Telescope at Paranal Observatory (VLT), and with the 
Low Resolution Imaging Spectrograph (LRIS) at the W. 
M. Keck Observatory (WMKO). Observational configu- 
rations and objects observed are listed in Table [2j 

Since at the time of observation, many of these sources 
were not classified, we often initially obtained only suf- 
ficient S/N to identify the redshift of an FSRQ, or to 
firmly establish the source as a BLL. Also, with the va- 
riety of telescope configurations and varying observing 
conditions, the quality of the spectra are not uniform: 
resolutions vary from 4 to 17 A, exposure times from 
180 s to 2400 s, and telescope diameters from 3.58 m to 
10 m. S/N per resolution element varies from 10 to > 300. 
In a number of cases follow-on observations with higher 
S/N and/or higher spectral resolution allowed us to more 
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carefully study confirmed BLL lacking redshifts. Here we 
discuss the most constraining spectrum or spectrum av- 
erage for each source, referring to this as the 'primary' 
spectrum. 

All spectra are taken at the parallactic angle, except 
for LRIS spectra using the atmospheric dispersion cor- 
rector, where we observed in a north-south configuration. 
In a few cases, we rotated the slit angle to minimize con- 
tamination from a nearby star. At least two exposures 
are taken of every target for cosmic ray cleaning. Typical 
exposure times are 2x900 s. 

2.3. Data Reduction Pipeline 

D ata reduction was perfor med with the IRAF pack- 
age (jTodvll 19861 : IValdes 1986) using standard techniques. 
Data was overscan (where applicable) and bias sub- 
tracted. Dome flats were taken at the beginning of every 
night, the spectral response was removed, and all data 
frames were flat-fielded. Wavelength calibration em- 
ployed arc lamp spectra and was confirmed with checks of 
night s ky lines. We employed an optimal extraction algo- 
rithm (Valdes 1992) to maximize the final signal to noise. 
For HET spectra, care was taken to use sky windows very 
near the longslit target position so as to minimize spec- 
troscopic residuals caused by fringing in the red, whose 
removal is precluded by the rapidly varying HET pupil. 
Spectra were visually cleaned of residual cosmic ray con- 
tamination affecting only individual exposures. 

We performed sp ectrophotom etric calibration u sing 
standard stars from lOkd ([1990ft and iBohlinl ([2007ft . In 
most cases standard exposures were available from the 
data night. On the queue-scheduled HET, and during 
our queue-scheduled VLT observations, standards from 
subsequent nights were sometimes used. At all other tele- 
scopes, multiple standard stars were observed per night 
under varying atmospheric conditions and different air- 
masses. The sensitivity function was interpolated be- 
tween standard star observations when the solution was 
found to vary significantly with time. 

For blue objects, broad-coverage spectrographs can 
suffer significant second order contamination. In par- 
ticular, the standard HET configuration using a Schott 
GG385 long-pass filter permitted second-order effects 
redward of 7700 A. The effect on object spectra were 
small, but for blue WD spectrophotometric standards, 
second order corrections were needed for accurate deter- 
mination of the sensitivity function. This corre ction term 
was constructed following iSzokolv et al.l ([2004) . In addi- 
tion, since BLL spectra are generally simple power laws, 
we used our objects to monitor second order contami- 
nation and residual errors in the sensitivity function. In 
extreme cases, we fit an average deviation from power law 
across all objects in a night, and treated it as a correction 
to our spectrophotometric calibrations. This resulted in 
excellent, stable response functions for the major data 
sets. 

Spectra were corrected for atm ospheric extinct i on us - 
ing standard values. We followed Ikrisciunas et al.l {1987) 
for WMKO LRIS spectra, and used the mean KPNO ex- 
tinction table from IRAF for P200 DBSP spectra. Our 
NTT, VLT, SOAR, and HET spectra do not extend into 
the UV and so suffer only minor atmospheric extinc- 
tion. These spectra were also corrected using the KPNO 
extinction tables. We removed Galactic extinction us- 



ing IRAF's de-r eddening function and the Schlegel maps 
(Schl egel et al.|[T998ft . We made no attempt to remove 
intrinsic reddening (i.e.: from the host galaxy). 

Telluric templates were generated from the standard 
star observations in each night, with separate templates 
for the oxygen and water line complexes. We cor- 
rected separately for the telluric absorptions of these two 
species. We found that most telluric features divided out 
well, with significant residuals only apparent in spectra 
with high S/N. On the HET spectra, residual second 
order contamination prevented complete removal of the 
strong water band red-ward of 9000 A. 

When we had multiple epochs of these final cleaned, 
flux- calibrated spectra with the same instrumental con- 
figuration, we checked for strong continuum variation. 
Spectra with comparable fluxes were then combined into 
a single best spectrum, with individual epochs weighted 
by S/N. 

Due to variable slit losses and changing conditions be- 
tween object and standard star exposures, we estimated 
that the accuracy of our a bsolute spectrophotometry is 
~ 30% (Heal ev~ et al.l l20Q8h . although the relative spec- 
trophotometry is considerably better. 

Fig. Set 1. Spectra 

2.4. General Trends 

To illustrate the principal trends in the BLL spectra 
we refer the reader to Figure [U By definition, the dom- 
inant component is a power-law. J0516— 6207, however, 
shows that after removal of a power-law weak, but broad 
C IV, C III, and Mg II features may occasionally be seen 
in high S/N spectra. Here the equivalent widths (< 1 
A) are sufficiently small to secure the identification as 
a BLL. However, should the continuum fade by ~ 10 x, 
this would be classified (at that epoch) as an FSRQ. Such 
transition objects support the idea of a Blazar contin- 
uum, rather than two distinc t pop ulations (Fos sati et al.l 
119981 : iGhisellini Tavecchiol 12008). In S12, we reported 
significant broad line measurements for 5 of our BLL, 
including J0516-6207. 

While the power-laws of most BLL are very blue, like 
J0516-6207, a few like J1849+2748, appear intrinsically 
flat or red, even after correcting for Galactic extinction. 
This may plausibly be a sign of a synchrotron compo- 
nent peaking near the optical, but might also indicate 
incomplete extinction correct ion, with residual r edden- 
ing caused by dust not in the lSchlegel et al.l (j 1998ft maps. 
It could also be intrinsic host extinction. 

The Galactic reddening can be very severe. For 
J2001+4352 (upper right) we show both the highly 
extincted, pre-correction spectrum and the blue post- 
correction power law. This source is in a direction of 
known high Ay — 1.75. Such extincted power-law spec- 
tra provide an excellent opportunity for ISM studies: The 
features seen after de-extinction and division by the best- 
fit power law (lower panel) are all interstellar in origin 
- Galactic H and K, Na I 5892, and a series of diff use 
interstellar bands, as described in lYuan fc Liul (j2012ft . 

For J0124— 0625 (upper left) the residual absorption 
features are intergalactic in origin. Redward of 3900 A we 
detect a number of metal-line systems, blueward one sees 
the onset of strong Lyman-a forest absorptions. These 
features determine a redshift z = 2.117 (one of the high- 
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Table 2 

Observing Configurations 



Telescope 


Instrument 


Resolution 
A 


Slit Width 
Arcseconds 


Objects 


Filter 


A 


A 


HET 


LRS 


15 


2 


41 


GG385 


4150 


10500 


HET 


LRS 


8 


1 


8 


GG385 


4150 


10500 


NTT 


EFOSC2 


16 


1 


31 




3400 


7400 


NTT 


EMMI 


12 


1 


1 




4000 


9300 


Palomar 200" 


DBSP 


5/15 


1 


4 




3100 


8100 


Palomar 200" 


DBSP 


5/15 


1.5 


5 




3100 


8100 


Palomar 200" 


DBSP 


5/9 


1.5 


42 




3100 


8100 


SOAR 


GHTS 


6 


0.84 


2 




3200 


7200 


VLT 


FORS2 


11 


1 


14 




3400 


9600 


VLT 


FORS2 


17 


1.6 


16 




3400 


9600 


WMKO 


LRIS 


4/7 


1 


90 




3100 


10500 


WMKO 


LRIS 


4/9 


1 


40 




3100 


10500 



Note. — For DBSP and LRIS the blue and red channels are split by a dichroic at 
5600 A; the listed resolutions are for blue and red side, respectively. 



est in our BLL sample). The lack of similar Lya forest 
absorption in many of our other high S/N, high resolu- 
tion spectra allows us to place statistical upper limits on 
the redshift as described in §3.31 

Finally at lower right we see two BLL with significant 
flux from the host. In J2042+2426, the galaxy provides 
about a third of the total flux and is easily visible in the 
raw spectrum. This is still safely a BLL, and we can 
measure the continuum contribution by a 'Non-Thermal 
Dominance' (see g5j here NTD=1.38). For J2055-0021, 
the host is swamped by the core synchrotron emission 
(NTD=47.5) and the galaxy features are visible only af- 
ter subtracting the best-fit power law, as in the lower 
right panel. The flux increase in the blue appears due to 
residual few-% fluxing errors (here, likely incomplete cor- 
rection for atmospheric extinction), rather than intrin- 
sic emission. Despite the careful calibration, such resid- 
ual fluxing issues persist in several spectra. However, 
the high-pass filtering described in §3.41 ensures that our 
measurements of, and bounds on, host galaxy flux are 
almost completely immune to such residual calibration 
artifacts. We find that a number of BLL show visible 
host galaxy componen ts, all consistent with giant ellipti- 
cals ( lUrrv et al.l l2QQQ). We discuss the flux distribution 
of these host galaxies in §3.71 

2.5. Individual Objects 

A number of BLL reductions required special treat- 
ment. For example a few objects clearly required changes 
to the Schlegel map Ay, so that the de-extinction re- 
sulted in clean power laws. For J0007+4712, Ay was 
increased from 0.3 to 0.8 and for J1941— 6211 from 0.3 
to 1.0. Conversely we decreased the Ay of J1603— 4904 
from 7.8 to 5.0 and J2025+3343, from 6.15 to 5.0. We 
checked th e recent recalibration of the Galactic extinc- 
tion maps (Schla ffv fc Finkbeinerll20lTt ). but did not find 
large changes, so these extinction features affecting our 
BLL are probably on scales below the map resolution. 

J1330+7001 was observed off of the parallactic angle — 
the ensuing drop in blue flux is not intrinsic to the sys- 
tem, and our power law is fit only redward of 5000 A. 
We thus remove the broadband residual in our power 
law divided spectrum in Figure [TJ J1829+2729 and a 
nearby star were spatially unresolved in our data. The 
presented spectrum is a composite of starlight and quasar 



light, which the significant emission features all identified 
as z = stellar or ISM features. 

In a few cases, objects previously cataloged as BLLs do 
not have sufficient S/N in our spectra for a definitive BLL 
classification. For J080 1+4401, we find that undetected 
broad lines could have an equivalent width as large as 9.5 
A; for J0209-5229 the limit is 5.6 A, for J1311+0035 the 
limit is 8.0 A and for J1530+5736 we could have missed 
lines as strong as EW=5.5 A. As higher S/N spectroscopy 
would likely confirm the archival BLL designations, we 
consider them BLLs for the purposes of this paper. 

Five of the BLL described here had high significance 
broad line detections and have already been described in 
S12; we re- measure these spectra here for a uniform BLL 
treatment. In J0847-2337 and J0430-2507, the flux and 
spectral index measurements differ from the S12 values. 
This is because in the present analysis we first subtract 
the host galaxy flux, and calculate the flux and spectral 
index of the remaining non-thermal component. In S12, 
no such correction was attempted. 

3. MEASURING BLL REDSHIFTS 

The opportunity to advance our understanding of BLL 
evolution with the large, flux limited Fermi sample is im- 
portant ([Ajello et al.ll2012bl ). Yet, despite the substan- 
tial telescope resources and careful analysis summarized 
above, many BLL did not yield direct spectroscopic red- 
shifts, due to the extrem e weakness of their emission lines 
(Sbarufa tti et al.ll2005bl ) and lack of clear host features. 
Therefore we collect here both the direct redshift mea- 
surements and quantitative constraints on the allowed 
redshift range for our observed BLL. 

3.1. Emission Line Redshifts 

We visually inspected all spectra for AGN emission 
line features, and host galaxy absorptions. Spectroscopic 
redshifts are meas ured by cross-correlat ion analysis using 
the rvsao package ([Kurtz fc Mink|[T998l ). We require one 
emission line to be present at the > 5<r level, and a second 
line present at the > 3a level — significances are measured 
by fitting a Gaussian template in the splot tool; we allow 
the width and amplitude of the Gaussian to vary, but fix 
the center at the redshift derived by rvsao's xcsao rou- 
tine. For this study, we limited our search to typically 
strong emission lines known to be present in some BL 
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Figure 1. Spectra of Fermi BLLs. Each object is presented twice - directly in the upper panel, and then with the best-fit power law 
'removed' (generally by division, but for composite spectra by subtraction). This residual is plotted in the lower panel. SDSS spectra, 
while discussed in Table [3] are not replotted in this Figure. The six sample spectra here illustrate major trends. Similar figures for the full 
BLL sample are available in the electronic edition of this journal. There we mark only the two lines used to secure z S pec] here other lines 
of interest are marked and the two qualifying lines are marked in red. z limits given are the most constraining limits presented — 2 digits 
are given for limits from non-detection of host galaxies; 3 digits, for intervening absorption systems. The marked absorption system is the 
best spectroscopic limit; in some cases, a stronger host galaxy limit is presented. See Table [3] for the precise spectroscopic limit. 



Lacs: Broad emission from C IV (1549, 1551 A), C III 
(1909 A), Mg II (2796, 2799, 2804 A), H 7 (4340 A), Rf3 
(4861 A), and Ha (6563 A) and narrow emission from 
[O II] (3727, 3729 A), [O III] (4959, 5007 A), and [N II] 
(6549, 6583 A). While other species exist in our spec- 
tra, these here listed are sufficient for definite redshift 
IDs. Velocities are not corrected to helio-centric or LSR 
frames. 

In many cases, spectroscopic redshift s are further de- 



termined by a significant (> 3 a) de tection of a host 
galaxy, as will be described in £13.41 A few redshifts 
require further note: In J0124-0625 and J1451+5201, 
we identify the redshift by a Lya and C IV absorption 
system at the onset of the Lya forest. In J0434-2015, 
we identify a single strong feature with [O II], consis- 
tent with weak Mg II and Ca H/K absorptions. For 
J1728+1215, we find strong Mg II, confirmed by [O II] 
at the same z in an archival spectrum. In J2152+1734, 
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we identify a strong feature with Mg II confirmed by a 
significant [O II] detection in archival spectroscopy. 

For a few objects only a single emission line was mea- 
sured with high S/N. In general we use the lack of oth- 
erwise expected features to identify the species and the 
redshift with high confidence. Nevertheless, a few red- 
shifts have some systematic uncertainty and are marked 
by a in Table [31 We briefly discuss these cases here. 
For J0007+4712, we derive a redshift from the clear on- 
set of the Lyman-a forest and report only two significant 
figures. In J0212+2244, we determine a tentative z from 
weak Ca H, K and g-band absorptions. For J0439-4522, 
we identify the one strong emission feature as C IV; in- 
tervening absorption excludes a Mg II identification, but 
a less likely C III identification at z ~ 1.45 is not conclu- 
sively ruled out. J0629— 1959 presents broad but weak 
emission at the redshift of the highest z metal line ab- 
sorption system (1.724). We thus identify this, tenta- 
tively, as the object's true z. For J0709— 0255, we iden- 
tify the strong feature at 9200 A with [O III] by the line 
shape; an [O II] identification at z ~ 0.84 is not ex- 
cluded. For J0825+0309, we find significant [O III] emis- 
sion at 5007 A (and possible, but not significant emission 
at 4959 A), at a z con s istent with an Mg II feature iden- 
tified in Stic kel et al.l (j!993f ). We find weak features in 
J1231+2847 at the SDSS z, but they have low signifi- 
cance. For J1312— 2156, we find a plausible Mg II fea- 
ture; this single line identification is in a small allowed 
redshift range (z ~ 1.6), other identifications for this 
line are spectroscopically excluded. In J 1754-6423, we 
tentatively identify emission at ~ 6300 A with Mg II — 
higher z redshifts are excluded by the lack of Lya forest. 
In J2116+3339's spectrum, a significant broad emission 
feature is identified with C IV, consistent with a weak 
bump in the far blue at Lya. Nevertheless a lower z 
redshift is possible if the purported Lya line is not real. 
J2208+6519 presents one strong, broad emission feature, 
tentatively identified as Mg II — a C IV identification at 
z ~ 1.8 is not excluded. 

3.2. Intervening Absorbers 

For some BLLs, the core light passes near an interven- 
ing galaxy on its way to Earth. At small radii one can 
encounter low excitation clouds in the galaxy's halo, giv- 
ing absorption doublets from Mg II at (2795.5, 2802.7) 
A. Larger impact parameters can sample C IV at (1548.2, 
1550.77) A. In some low excitation (i.e.: Mg II) systems, 
we also see absorption from Fe II at (2344.2, 2374.4, 
2382.7, 2586.6, 2600.1) A. Finally, for our highest red- 
shift BLL we Lyman-a absorption at 1215 A for the 
metal line systems, as well as onset of the Lyman-a for- 
est. 

For unsaturated absorptions, the doublet ratio for 
Mg II and C IV is 2:1, with the blue li ne dominant. 
In saturated absorptions, the r atio is 1:1 ([Nestor et al.l 
120051 : iMichalitsianos et atlll988h . 

We visu ally search all spect ra for candidate doublets, 
and follow [Nestor et ahl ([2005) in employing a quantita- 
tive test of the significance of each candidate. We used 
a two Gaussian template with wavelength spacing scal- 
ing with z, but free amplitudes, and fit the equivalent 
width and error of each component, using the splot tool 
in iraf. For the candidate to qualify as a detection we 



require the stronger (bluer) line to have > 5<r signifi- 
cance, and the second line to have > 3<r significance. In 
a few cases, one component of an otherwise strong dou- 
blet was affected by skylines, telluric features or cosmic 
rays. In these cases, another expected feature from the 
absorption complex (e.g.: a Fe II line) detected at > 3a 
qualified the system. We further require the doublet ra- 
tio to be consistent (within errors) to a value between 
2:1 and 1:1. 

Our principal goal is not an absorption line study. 
Thus we concentrated on the longest wavelength (high- 
est z) candidate system and measured sufficient lines to 
confirm the z (i.e: two significant lines). After validation 
we continued to search for higher z until no candidates 
passed the significance test. We therefore believe that 
we have found the highest z absorption system in each 
of our spectra strong enough to reach the ba/3a criteria 
above. 

Since we see the onset of the Lyman-a forest in our 
highest redshift objects, the red end of the forest provides 
a strict lower limit on redshift. This can be higher than 
that inferred from the reddest metal line system. 

We list these spectroscopic minimum z's as z m { n , when 
available, in Table 02 

3.3. Redshift Upper Limits 

We can use the absence of Lyman-a absorptions to 
provide statistically-based upper limits on z for all BLLs 
without redshift. The exclusion z max depends on the 
spectral range, resolution and S/N of the particular ob- 
servation, but is generally 1.65 < z < 3.0. 

To quantify the upper boun d, we need the expe cted 
density of Lya forest absorbers. iPenton et al.l ([2004) find 
that for rest EW > 0.24 A, dN/dz ~ 40 at z = 1.6, 
varying wi th redshift as log dN/ dz oc 1.85 log (1 + z). 
We follow IWevmann et al.l (|1998[ ) for the EW scaling: 

dN/dz (X e -(EW test -0.24)/0.267 m Ag the g/ N and resolu _ 

tion of our data vary, we generally measured a conser- 
vative uniform 3a sensitivity for narrow Lya absorptions 
100 A from the blue end of our spectra. Typical EW lim- 
its in this range were ~ 0.2 — 1.0 A. Given this density 
we solve for the Az range giving a Poisson probability of 
0.32 (ie: la) for detecting no absorbers, obtaining: 

/ \ . \-i.85 

AZ = 0.167 • e (^est-0.24)/0.267) . / Vn \ (1) 

where A m i n is the effective blue end of the spectrum (gen- 
erally 3150 - 4200 A), and EW is the measured equivalent 
width limit. Thus, we infer a maximum source redshift 
z max = (A m in/1215) — 1 + Az. In some cases, the S/N 
is too low at the blue end of the spectrum. We then 
measure EW limits closer to the sensitivity peak of the 
spectrum. Of course, with a larger A m i n for the effective 
spectrum end Equation [1] gives less constraining upper 
limits. 

If the actual blazar redshift z is very close to £ max as 
estimated above, then its UV radiation may photo- ionize 
Lya clouds along the line of sight, postponing the onset 
of the forest and artificially lowering z max . In practice, 
the effect is usually small (Az 

max ^ 

0.01 - 0.02) except 
for large z max wh en our bound is ge nerally not of in- 
terest. We follow iBajtlik et al.l (1988) to estimate this 



7 



'proximity effect' correction, by computing 



4ttJ u (1 + z) 



(1 + ^) 



1/2 



-i 2 



(1 + 2)1/2 _ (1 + ^ ya )l/2 

(2) 

where f v is the blazar flux at the absorption Lyman 
limit (at ^Lya) and J v = 10 -21,5 ergcm -2 s _1 Hz -1 sr _1 
is the cosmic ionizing flux, conservatively estimated for 
our redshift range. We compute f v from the power law 
fit in Table [3j and increase the blazar redshift z until 
(jo(z\, YOi ) < 1. We thus adopt these revised z max = z and 
quote these corrected upper limits in Table 02 

3.4. Host Galaxy Fitting 

It has been claimed that BL Lac objects are hosted by 
giant elliptical galaxies with bright ab solute magnitude 
- M R = -22.9 ± 0.5 in our cosmology (jUrrv et all 120001 : 
iSbarufatti et al.ll2005bf ). 

If we adopt the common assumption that these are 
standard candles, we can estimate the redshift of the BLL 
by detecting such galaxies. In imaging studies, one looks 
in the wings of the BLL for the host galaxy flux, and 
compa res that to the standard candle flux at various red - 
shifts (jSbarufatti et alJl2QQ5at iMeisner fc Romanill2010f ). 
In spectroscopic studies, one typically looks for individ- 
ual absorption features (i.e.: H, K, g-band). One can 
also use the lack of such lines as evidence that the BLL 
is at higher redshift (|Sbarufatti et al. 2005b; S haw et aTl 
[20(1). 

With high S/N spectra, however, one can obtain more 
stringent limits by using the entir e elliptical template 
rather than just one (or a few) lines. iPlotkin et al.l ([2010) 
developed a technique of fitting for host galaxies in SDSS 
BLL spectra. We expand here on that method for our 
more heterogeneous spectra. 

Our spectra come from a variety of spectrographs in 
disparate observing conditions and we find low frequency 
systematic fluctuations in many of the fits. These are 
likely caused by imperfect spectrophotometric flu xing 
and second order contamination as discussed in £12.31 
These effects can dominate over real galaxy features. 
Using SciPy's Signal Processing routines 4 , we construct 
a bandpass Kaiser window from 220 A to 1.5 x the 
Nyquist frequency. We apply that window as an effec- 
tive high pass filter both to our spectra and to the tem- 
plates, to mitigate this l ow- frequency noise before fitting 
(Kaiser fe Sc hafer 1980h. 

We test possible redshifts Z{ on a grid scaled to the 
spectrograph resolution with constant spacing in logz. 
This grid is thus z { = (^ + l) i - 1, where AA is 
the pixel scale of the spectrograph. For each trial Z{ 
we fit the power law F v oc v a index and flux and 
the amplitude of a redshifted elliptical template. This 
host template is generated from the PEGASE model 
(|Fioc fc Rocca-Volmer ange 1997) tables and evolved to 
low z from z = 2, as in lO'Dowd fc Urrvl([2005l ). For uni- 
formity, we here use the same template for all Zi, and do 
not perform any evolution corrections. 

Our fit minimizes y 2 with three free variables at each 
trial Zi. We employ the scipy.optimize.leastsq routine 

4 Documentation and more in formation available at 
http: / /docs, scipy.org/doc/scipy/reference/signal. html 



based on a Levenberg-Marquardt fitter. 

To model host slit losses, we assume an r = 
lOkpc de Vaucouleurs profile with Sersic index 1/4 
([O'Dowd fc Urrvl l2QQ5h and account for the individual 
observations' slit width and seeing profile (measured 
from the core full width at half maximum, FWHM). 
Since we have employed optimal extractions of the BLL 
spectra, our effective aperture along the slit varies, but 
we estimate a typical width of ~ 2x the spatial FHWM 
achieved during our spectral integration. Accordingly we 
estimate host slit losses through a rectangular aperture 
of the slit width x twice the spectrum FWHM. Inferred 
host fluxes are corrected for these slit losses. 

Results of sample fits are shown as the blue dots in Fig- 
ure [2] where the fit amplitude of the host galaxy template 
is plotted against trial redshift. 

3.5. Power Law Fit 

We report the power law fluxes and spectral indices of 
the best fit to the de-extincted spectrum in Table [3l The 
flux is given in units of Log 10 _28 erg cm _2 s _1 Hz _1 as ob- 
served at 10 14 T Hz (~ 5980 A), the center of our typical 
spectral range. The index a is measured F v oc v a . These 
values may be combined with multi-wavelength data to 
study the continuum SED of the blazars in our sample. 
Since the statistical errors on the fit are, in general, un- 
physically small, we follow S12 in estimating errors on 
the spectral index by independently fitting the red and 
blue halves of the spectrum. Note that large errors bars 
generally indicate a relatively poor fit to a simple power 
law rather than large statistical errors. The statistical 
errors on the F v amplitude are also small; we convolve 
these with ou r estimated 30% overall fluxing uncertainty 
(Hea lev et al.ll2008h . which dominates in nearly all cases. 

For objects with high significance (> 3<r) detections 
of galaxies, as described in §3.7| we report the best fit 
power law from the simultaneous power-law/host fit in 



When we have observed objects at multiple epochs, we 
also fit a power law to the other, non-primary spectra. 
These fluxes vary substantially, some by more than 10 x. 
In Figure [3j we show the distribution of / m ax//min flux 
ratios. This is well-described by an a = 1.29 power law. 
Epochs from our fiducial spectra are listed in Table 3. 

3.6. Testing the Standard Candle Assumption 

BLL with a redshift and a secure (> 3a) host de- 
tection can be used to test the uniformity of the host 
luminosities. There are 59 such BLL in our sample. 
We derive synthetic i?-band magnit udes by applying a 
Kron-Cousins R filter to our spectra (jMeisner fc Romanil 
l2010h . The results are shown as a histogram in Figure 
H We find < M R >= - 22.5, down ~ 0.4 magni tudes 
from M R = -22.9 found in ISbarufatti et al.l (|2QQ5b[ ). We 
find a similar spread in luminosity (~ ±0.5 magnitudes). 
When we separate the host measurements for lower-peak 
sources (LBL+IBL) we find that they have a median 
luminosity ~ 0.3 magnitudes fainter than that of our 
HBL. Unfortunately, we do not have sufficient LBL+IBL 
hosts to test for su ch differences at high significance. 
Past studies differ: lUrrv et all J20OO) found no signif- 
ica nt offset in the host m agnitudes of HBL and LBL, 
but ISbarufatti et al.l (|2005bf ) noted that higher peak HBL 
tend to have more luminous hosts. 
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Figure 2. Best fit host galaxy fluxes as a function of trial zi (blue dots). Our estimate of the 2a local systematic flux error is shown 
by the blue line (see text). The red bands give the la range about the expected host flux for Mr = —22.5 db 0.5. Vertical green dashed 
lines give our z m i n (— 22.5) limit on the redshift. For comparison, solid green lines give the spectroscopic redshifts, when measured. For 
J2055— 0021 and J0148+0129 these are consistent; the former has a high significance host detection (Figure while for the latter the 
higher z redshift is from emission lines. J1534+3715 is one of 9 objects with limits inconsistent with spectroscopic redshifts. Here, the host 
galaxy is sub-luminous (Mr = —21.87 db 0.16) and thus (just) missed by this technique (Prob = 0.15). 



Two of the high significanc e host galaxies have imag - 
ing magnitudes reported in Sb arufatti et al.l ([20 05b). 
For J1442+1200 we m easure M R = -22.99 ± 0.14; 
iSbarufatti et al.l (j2QQ5bh reported Mr = -22.77. For 
J1428+4240, we meas ure M R = -22.78 ± 0.13; 
ISbarufatti et al.l (j2005bf ) find M R = -22.75. These are 
consistent within measurement errors, a good check of 
our slit-loss corrections and magnitude estimates. 

Overall, the LAT BLL sample thus represents a fainter 
host population of than those studied in previous work. 
Conceivably a higher (LBL+IBL) /HBL ratio in our sam- 
ple causes part of the difference (although we remain 
HBL dominated). However, we suspect that our rather 
exhaustive 8-m class campaign, skipping most objects 
with redshifts already in the literature, selects for fainter 
host galaxies than in the past. Thus we may be probing 
fainter on the true host luminosity distribution; the BLL 
for which we were not able to provide host detections may 
then be similarly under-luminous compared to previous 
studies. A true evaluation of the intrinsic host luminosity 
distribution, as well as any dependence on subclass type, 
will require a careful assessment of the parent population 
(7-ray) and host detection selection biases. 

In the rest of this section, we conservatively adopt our 
Mr = — 22.5 ± 0.5 estimate. We do also report (Table 
[3]) more aggressive redshift limits based on the common 



assumption Mr = —22.9 ± 0.5, for more direct compari- 
son to previous work, but we recommend adoption of the 
less stringent Mr — —22.5 ± 0.5 redshift constraints. 

3.7. Lower Limits from Non- detections of Host Galaxy 

We use the results of the fitting in §3.4l and the calibra- 
tion in §3.6l to constrain the redshift of the host. At each 
trial redshift, our fitter reports a flux (/ ± af) for the 
host galaxy. This is to be compared with the model flux 
from the redshifted standard candle elliptical template 
(/m±A/ m ). 

While we have greatly decreased the effect of the low 
frequency noise in our spectral fits using the high-pass 
filter, we still find that the statistical errors on the fit 
host galaxy fluxes are unrealistically small; these flux 
estimates remain dominated by systematic effects. We 
therefore construct an error estimate Af for each Zi by 
measuring the dispersion of flux estimates for nearby red- 
shift bins. This is computed from a sample of the 30 
nearest f(z), skipping 5 bins on each side of our test Z{ 
to minimize high pass correlation. After a-clipping, the 
fit flux distributions are well behaved and we use these to 
compute a 2-a upper limit (scaled from the rms) at each 
Zi. Vectors of these upper limits are shown by the jagged 
blue lines in Figure [21 These have captured the local ef- 
fective noise quite well and so we adopt these vectors as 
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Figure 3. Histogram of maximum flux ratios. For each object 
where we have collected multiple spectra, we compute the flux ratio 
between the power law components of the brightest and faintest 
observed epoch. These are plotted as a log-log histogram. The 
dashed line is the best-fit power law with index a = 1.29. The inset 
histogram shows the distribution of time between observations - 
typically ~ 1 year 
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Figure 4. Measured BLL host absolute magnitudes (R equivalent 
at z = 0). Black histogram: all BLL hosts with spectroscopic red- 
shift. A weighted x 2 fit gives Mr = —22.5 ±0.5, our best estimate 
for the luminosity of BLL host galaxies. Red: hosts with high sig- 
nificance (> 3cr) detections. Green and Blue sub-histograms show 
high significance LBL+IBL and HBL hosts, respectively. 

effective 2<r confidence limits. 

Fit fluxes well above these 2a limits denote likely host 
detections. Indeed, we found that this automatic pro- 
cessing was quite effective at flagging candidate Z{ for 
host detection. Here, however, we focus on how well a fit 
flux / with local effective error Af can be used to exclude 
a host of the expected magnitude /m at the test redshift 



z. Thus we compute a probability that the fit flux / and 
error Af are consistent with the expected model flux Jm 
and its uncertainty at the given z as 

/>oo 

Prob(z; f, Af) = C G(f' M , Af' M )-G(f, Af)df' M (3) 



where G(x,a x ) is a Gaussian of width a x centered at x. 
The normalization C is chosen such that Prob = 1 for 
Im = (i.e.: any / is acceptable for a model of zero 
flux) . This is a conservative choice as it does not exclude 
over-luminous hosts. For example, when we assume a 
model Mr = — 22.5 ± 0.5 this probability also finds any 
/ consistent with Mr = —22.9 ± 0.5 to be acceptable. 
When the fitter returns an unphysical negative /, we 
evaluate the probability for / = and the local Af. 

This probability becomes substantial for z near a good 
candidate redshift. It also grows as the sensitivity of 
our host search drops at large z. We thus calculate a 
minimum redshift (z rnin ^_ 2 2.b)) as the lowest redshift for 
which Prob > 0.17 (i.e.: la). We list these values in Ta- 
ble[3l For comparison, we also give z m i n (-22.9) calculated 
in the same fashion, assuming a model Mr = — 22.9±0.5. 
Comparison between the spectroscopic detections and 
z m [ n suggest that the Mr = —22.5 value is most con- 
sistent with observed detections, and lower bounds (as 
expected from §3.6). We recommend use of these conser- 
vative lower limits. Note also that the vector Prob(z), 
once normalized with an appropriate prior and truncated 
at z max , can be used as a PDF for the BLL redshift. 

3.8. Redshift Distribution 

As seen in Figure [5j archival redshift measurements for 
BLL are dominated by low values (z = 0.23). Our new 
spectroscopic redshifts have extended the population to 
higher z, with some objects' redshifts at z ^> 1. Still, the 
objects with redshifts remain dominated by low z. In the 
new spectroscopic redshifts we find z = 0.33. We believe 
there is a significant bias to low redshift in both of these 
samples, as the weak low EW emission or absorption fea- 
tures of our typical BLLs with known redshift are easier 
to detect at low z. 

In Figure [5] we also show two sets of redshift lower lim- 
its. For every object in our sample, we can derive a host- 
detection limit (z = 0.41). These lower limits on redshift 
are still biased low: as evident from Figure [2] we are most 
sensitive to galaxies at low z. Nevertheless, they suggest 
that these objects are not consistent with the spectro- 
scopic redshifts (the K-S test gives probability < 10 -11 
of consistency with archival redshifts). The absorption 
line limits we have for some objects (described in §3.2p 
give further evidence for a population of BLLs at high 
redshift (z > 1). Together, these results strongly im- 
ply that previous BLL studies suffered important biases 
due to shallow samples with large redshift incomplete- 
ness preventing detection of bright, but high z BLL. 

The inset shows the spectroscopic redshifts and the 
redshifts limits for high-peaked (HBL) and lower-peaked 
(LBL+IBL) sources classified in 2LAC. We see that the 
lower-peaked detections extend to higher z, as might be 
expected if these sources are more luminous and have 
a less dominant synchrotron continuum. However, the 
open histograms of limits remind us that both sub-classes 
still suffer substantial redshift incompleteness, and the 
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Figure 5. BLL redshift and redshift lower limit distributions. 
Within each 5z = 0.2 bin, we show (left to right) the Archival z, 
our new spectroscopic z, lower limit z's from host fitting and lower 
limit z's from intervening absorption line systems. The limits show 
that the spectroscopic redshift samples, particularly the archival 
sample, are selection biased to low z. The subframe shows the red- 
shifts (filled histograms) and lower limits (open histograms) for the 
LBL+IBL (green) and HBL (blue). Again, both limit histograms 
extend to higher z. The LBL+IBL sample extends somewhat fur- 
ther than the HBL. 

missing redshifts for both sub-classes extend substan- 
tially higher than those in hand. A re-appraisal of the 
BLL population, properly including the new redshifts, z 
constraints and remaining selection biases is needed to 
test whether either subclass is still consistent with nega- 
tive cosmological evolution. 

4. BLACK HOLES AND HOST GALAXIES 

The masses of the central black holes provide impor- 
tant insight into the cosmic evolution of various AGN 
classes. These are most easily estimated by the virial 
technique (cf., iShen et all |2Q11[ ). In S12, we adopted 
this method to give mass estimates for the Fermi FS- 
RQs. For BLL, the lack of high S/N broad line mea- 
surements precludes such estimates. However, we have 
measured a number of host magnitudes in H3A\ since 
these are ellipticals, this is all 'bulge,' and we can apply 
an M — L relat i on to estimate the hole mass. We follow 
IGultekin et all (|2QQ9h : 



log( ^ 



S.95±0.11)+(l.ll±0.18)log 
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where M is the black hole mass, and Ly is the luminosity 
in a V filter [log(Ly/L ,y) = 0.4(4.83 - M° bulge )]. To 
convert our fit template spectrum amplitudes to consis- 
tent V magnitudes, we integrat e our template spe ctrum 
over the Hubble F555W filter (jLauer et al.ll2005l ) as in 

The masses from Equation |4] are plotted as circles 
in Figure [6] When the sub-class is known (HBL or 



LBL+IBL) we fill these in with blue or green, respec- 
tively. For comparison, we show the la spread of virial- 
esti mate BH mas s es from optically selected SDSS quasars 
from She n et ahl (j2011l ) (gray band) and masses of the 
Fermi FSRQs in S12 (red points). Of course, if BLL 
hosts really are standard candle ellipticals, then the 
M^uige — M. relation implies constant black hole masses. 
The masses corresponding to the standard Mr = —22.9 
and our revised —22.5 are shown by dashed lines. 

Interestingly, our BLL M. estimates increase with z 
much as the optical QSO or FSRQ. Of course, we only 
plot high significance host detections here, and low lumi- 
nosity hosts at high z are increasingly difficult to detect 
(unless the core luminosity decreases). Accordingly, as 
for the QSO, we suspect that the bulk of this trend is 
due to selection effects. In the case of the FSRQ, S12 
argued that the offset to smaller black hole mass was 
at least partly due to a preferentially polar view of an 
equatorially flattened broad line region, with the projec- 
tion decreasing the observed kinematic line width and 
the average virial mass estimate. Like 7-ray selected 
FSRQ, BLL are Dopple r-boosted along our line of sight 
(|Urrv fc PadovanH Il995f ) . However since the host flux 
is nearly isotropic, we expect little alignment bias in our 
M. estimates. Thus, it is unclear whether the BLL offset 
to larger black hole masses is real or selection dominated. 

In a study of BLL hosts detected in the SDSS 
iLeon-Tavares et aD (|2011l ) found no significant difference 
between the masses of the central black holes of HBL 
and LBL. In figure [6] the HBL masses are however biased 
upwards with respect to the lower-peak BLL black hole 
masses. This is of course just a restatement of the off- 
set in host luminosity seen in Figure HJ Unfortunately, 
we cannot claim that this is a physical difference as the 
trend is precisely what would expect from selection bias 
if HBL have brighter non-thermal cores. 

Ideally we could use these black hole masses to ex- 
plore the relationship between the BLLs and the general 
QSO population. Large black hole masses, if not induced 
solely by selection bias would imply a late stage of AGN 
evolution. The black hole mass is often compared to the 
source luminosity to characterize the state of the accre- 
tion in Eddington units. However, with the exception of 
the few BLL for which we see broad lines (which seldom 
have a significantly host detection), the observed flux is 
so dominated by beamed jet emission that quoting the 
accretion luminosity in Eddington units is not feasible. 

5. NON-THERMAL DOMINANCE 

In S12, we introduced the non-thermal dominance 
(NTD) as a quantitative measure of how much the op- 
tical is contaminated by non-thermal synchrotron emis- 
sion. We here extend that analysis to BLLs. 

For most BLL, the dominant 'thermal' contribution is 
the host galaxy, not the big blue bump. We therefore set 
NTD 

= -^core / -^host where both fluxes are measured at 
5100 A - the same wavelength as the H/3 continuum mea- 
surements for FSRQs. The wavelength choice is impor- 
tant for BLL NTD measurements, since the host galaxy 
is much redder than the continuum-dominated core; an 
NTD measurement just above (redward of) the 4000 A 
break would typically give values > 4x larger. Measure- 
ments below the break would, of course, diverge even 
more. 
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Figure 6. BLL black hole masses plotted as a function of redshift 
(circles; blue fill HBL, green fill LBL+IBL). The gray band shows 
the l a spread of the optically selected SDSS black holes ( Shen et al. 
2011). Red error bars give the la spread Fermi FSRQ black hole 
masses (S12). The estimated BLL black hole masses are higher 
than for the two other AGN types. The upper (lower) blue lines 
show the black hole mass expected for the standard candle bulge 
at M R = -22.9(-22.5). 

As noted BLL are highly variable, so that their NTD 
changes over time. Our BLL can vary by over an order of 
magnitude in the optical (S12, Fig. 3). In this study, our 
primary spectra for variable objects were generally drawn 
from the epochs giving the highest signal to noise on weak 
emission or absorption features. As this favors low core 
fluxes, our primary epoch biases the results towards low 
NTD. 

In Figure [7J we plot the NTD against 7-ray spectral 
index. The values for BLL with measured galaxy flux 
are plotted, along with lower limits for BLL with red- 
shift, but no significant host detection. The histogram 
at right shows the spectral index distribution of BLL 
with no redshift detection (and unknown NTD). We ex- 
pect these BLL to have higher NTD on average, as the 
larger core flux makes redshift determinations more chal- 
lenging. When known, we indicate whether the BLL are 
high-peaked sources (blue) or lower-peaked (green). The 
FSRQs from S12 are plotted as red triangles. 

The most striking trend in this figure is the vertical 
color separation. This is the well-known result that 7- 
ray spe ctral index hardens from FSRQ through LBLs to 
HBLs (|Ackermann et aD l2Qllh . We also see the defin- 
ing characteristic of the BLL, increased continuum dom- 
inance with respect to the FSRQ. However, the NTD 
trend does not appear to continue through the BLL: 
harder spectrum BLL do not in general show increased 
NTD. Indeed the highest NTD seem to be associated 
with LBL/IBL measurements and lower bounds. A plau- 
sible interpretation is that the LBL/IBL are more lumi- 
nous and hence visible to higher z, where we are less 
likely to detect a host galaxy. 

6. CONCLUSIONS 

We have dramatically increased the redshift complete- 
ness of this largest ever 7-ray selected BLL sample; 
~ 44% of these BLL now have spectroscopic redshifts, 
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Figure 7. The NTD plotted against 7-ray spectral index. The 
BLL are plotted as dots; solid if L 7 > 10 44 ; empty otherwise. 
Broad line objects from S12 (typically FSRQs) are plotted as red 
triangles - circles surround those in a BLL state at the primary 
spectrum epoch; squares, those known to be BLL at a different 
epoch. BLL with known z, but no thermal (host galaxy or broad 
emission line) detection are plotted as lower limits. BLL without 
known z are shown by the histogram at right. Blue symbols indi- 
cate HBL, green LBL+IBL. Except for the NTD increase for the 
BLL population relative to the FSRQ, no clear trend with spectral 
index is seen. 

and ~ 95% have at least a strong lower bound on the 
redshift. These z constraints show that the subset with 
actual spectroscopic redshifts is strongly biased to low 
z. Although we find that the measured redshifts for low- 
peaked BLL (LBL+IBL) do extend higher than those for 
the HBL with the highest synchrotron peak, our set of 
lower limits for both subclasses extend to yet higher z 
than the spectroscopically solved objects. Thus the ac- 
tual redshifts for all BLL are biased low compared to a 
flux-limited parent population of 7-ray BLL. This must 
be taken into account in any study of BLL evolution over 
cosmic time (Ajello et al, in prep). 

Many of our redshift limits rely on the common as- 
sumption that BLL hosts are standard candles. Our ef- 
fort to re-calibrate the standard luminosity for the Fermi 
BLL has resulted in fainter absolute magnitudes and, 
hence, more conservative minimum z. This implies that 
the standard candle assumption deserves further study, 
which could be best prosecuted by obtaining more high 
spatial resolution images of BLL with known redshifts. 
Our study provides a large increase in the spectroscopic 
redshifts, a useful precursor to such work. We find lim- 
ited evidence that HBL have more luminous hosts than 
LBL/IBL. Whether this is intrinsic or a selection effect 
in the presence of a brighter, harder continuum is not yet 
clear. 

Of course true spectroscopic redshifts are always 
preferable for uniform population studies. However, we 
suspect that much higher completeness will be difficult 
to attain, and will likely require novel observational tech- 
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niques, as significantly more time on 8 m class telescopes 
is both expensive and likely to provide only marginally 
better results. In the interim, the best hope for progress 
lies with careful correction for selection effects in the 
present sample. Because of the tight correlation of syn- 
chrotron peak frequency with LAT-measured spectral in- 
dex and because of the 7-ray index dependence of the 
LAT sensitivity, study of the relative population and evo- 
lution of different BLL subclasses will require correction 
for 7-ray selection effects as well as possible biases in the 
redshift determinations themselves. While we do not at- 
tempt such study here, our improved completeness will 
help in understanding these selection effects. 

We find that BLL black hole mass estimates (at a 
given redshift) are larger than those for optically selected 
quasars or Fermi FSRQs. Associated with the apparent 
host luminosity differences, the present detections sug- 
gests that HBL host the largest mass black holes. We 
cannot at present tell whether these trends are a true 
difference in the black hole populations or luminosity- 
driven selection effects. 

BLLs are, by definition, non-thermally dominated 
(NTD > 1). We find the BLL population to have signif- 
icant higher NTD than the Fermi FSRQs. This claim is 
conservative, as the objects without redshift are likely to 
be at even higher NTD. Since BLL have a range of 7-ray 
spectral index harder than those of FSRQ, it is perhaps 
surprising that this trend does not hold within the BLL 
class: NTD is not, on average, higher for the hardest- 
spectrum BLL. It seems likely that, in this respect at 
least, BLL are a distinct population from FSRQ, and 
not just the hardest spectrum objects on a continuum. 
In this case we might expect NTD to be controlled by the 
precise accident of the jet Doppler boosting. This accords 
well with the idea that NTD can vary for an individual 
source as the jet angle or effective Doppler factor vary, 
with relatively little change in the GeV spectrum. Multi- 
epoch, multiwavelength studies of the brighter BLL with 
a range of 7-ray hardness are needed to test these ideas. 

We thank the referee for a careful reading and com- 
ments that helped us add materially to the paper. 
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Table 3 

BLL Spectral Properties 
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RA 


Dec 
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z 
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—7.774111 
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1.55+0 


11 


— 1 399+0 003 








0.33 


0.50 


1.84 


BLL 


IBL 


P200 
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T00D7 8+4713 


1 QQQ862 


47.202135 


70007+4712 


1.33+0 


14 


— 1 431 +0 430 


2.1: 


S 
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11 
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11 


-0.776+0.010 








0.79 
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J0012.9-3954 


3.249804 
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Note. — Table[3]is published in its entirety in the electronic edition of this journal; A portion is shown here to show its form and content. 
Method for z ID — B, broad emission lines; N, narrow emission lines; G, host galaxy features, S, special case — see §3.1 
Spectroscopic lower limits (i.e.: From intervening absorption systems) 



